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Abstract. We present a quasar candidate identification 
technique based on multicolor photometry. The tradi- 
tional multi-dimensional method (2 x N dimensions, 
where N is the number of the color-color diagrams) is 
reduced to a one-dimensional technique, which consists 
in a standard fitting procedure, where the observed spec- 
tral energy distributions are compared to quasar simu- 
lated spectra and stellar templates. This new multicolor 
approach is firstly applied to simulated catalogues and 
its efficiency is examined in various redshift ranges, as a 
function of the filter combination and the available observ- 
ing time for spectroscopy. We conclude that this method 
is better suited than the usual multicolor selection tech- 



niques Ito quasar identification, especially tor high-rcdshitt 
quasars. The application of the method to real quasar sam- 
ples found in the literature results in an efficiency com- 
parable to the one obtained from the use of color-color 
diagrams. The major advantage of the new method is the 
estimation of the 'photometric redshift of quasar candi- 
dates, enabling, in almost all cases, spectroscopy to be 
targeted to best suited wavelength ranges. 
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of the sample, and to achieve a maximum completeness 
over a redshift range as large as possible. 

These surveys are essentially based on the photometric 
preselection of quasar candidates, in order to optimize the 
efficiency of telescope time. Preselection makes the com- 
pleteness questionable both at low and at high redshift 
(see discussion in e.g. Miller & Mitchell, 1998). Even mul- 
tiple color-color diagrams ( Hall et al., 1996b| ) still reveal 
a bias in the range 2.2 < z < 3, exactly where a revers- 
ing of evolution is thought to occur, making completeness 
specially important. At low redshift, it has been shown 
that the usual morphologic prerequisite for quasi-stellar 



candidates results in a noticeable deficiency (Wolf et al 



1998). This trend may increase the apparent evolution of 



quasars. 

The present paper describes a joint method both for 
distinguishing between quasars and stars/galaxies by their 
photometry and for obtaining an estimate of the photo- 
metric redshift of the quasar candidates as well. In we 
present the current state of the art in the optical search for 
quasars. §[| gives a brief description of our method, which 
has been first validated with simulations, as shown in §|J 
before being applied to real samples (§||). The last section 
of this paper (fQ) contains a summary of our conclusions 
and a brief discussion. 



1. Introduction 

Nowadays, very high redshift galaxies and quasars are be- 
ing found. However, quasars still offer the opportunity of 
a (relatively) easy construction of high-z samples, due to 
their high luminosity, quasi-stellar shape (allowing accu- 
rate photometry), and to their emission lines (allowing 
spectroscopical identification). Therefore, many large op- 
tical quasar surveys are in progress or in project. Even if 
high energy emission is quite a common quasar property, 
X-ray samples still need the detection of emission lines in 
the optical domain for the purposes of redshift determina- 
tion. To improve the efficiency of all quasar studies, one 
has to increase both the size and the upper redshift limit 
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2. State of the art in the optical search for 
quasars 

One of the most common problems in optical surveys is 
the selection criteria one has to apply in order to dis- 
tinguish between stars and quasars, since both categories 
of objects have (in most cases) a point -like morphology. 
Quasars are spatially resolved only at low redshifts and 
only when high quality imaging data are available. Selec- 
tion of QSO candidates among stellar objects from pho- 
tometric data is becoming a standard method. The UVX 
technique (UV excess) separates stars from quasars with 
a redshift z < 2.2. This technique is still currently being 
used, e.g. in the Chile-UK quasar survey, and the AAT 
2dF QSO Survey. 
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To avoid both the redshift limitation (z < 2.2) and the 
bias towards blue objects, other equivalent techniques are 
also being applied (for example BRX), with quite a high 
efficiency in other redshift ranges. Many Multicolor Sur- 
veys use an increased number of filters, especially towards 
the red, with special care paid to high-z quasars: examples 
are found in Boyle et al. (f99f), Kennefick et al. (f997), 
Osmer et al. (f998) and the references therein, Jarvis & 
Mac Alpine (1998) and the DPOSS (Digital POSS 11). 
In all these cases, the efficiency never exceeds 35% (~ 3 
candidates selected per one real quasar), depending on 
magnitude and redshift. 

For all the selection methods described above, the 
working space is the N x 2D-manifold of N color-color 
planes. The locus of quasars at different redshifts in each 
of these planes is compared to the locus expected for Main 
Sequence stars. For quasars, the allowed multi-color space 
is defined through simulated template spectra, set to dif- 
ferent redshifts. Basically all the selection procedures use 
the positions of candidates in the multi-color space ei- 
ther to compare with the expected locus for quasars or 
to compute the distance to the Main Sequence (Newberg 
& Yanny, 1997; Sloan Digital Sky Survey - Fan, 1999; 
Krisciunas et al., 1998). Quasar candidates are among 
the objects at the largest distances with respect to this 
Main Sequence "Snake" . However the selection of can- 
didates with quasar spectra is basically a ID problem: 
the fitting between the observed spectral energy distribu- 
tion (SED) and the equivalent one obtained from tem- 
plates (quasars/stars), all of them computed within the 
same photometric system. Two quasar surveys make use 
of such a fit to select candidates and to estimate their pho- 
tometric redshifts. The Calar Alto Deep Imaging Survey 
(CADIS) uses specialized filters and a multicolor classi- 
fication algorithm to distinguish between stars, galaxies 
and quasars, and the identification of high-z quasars in 
particular is based on the photometric redshift method 
(Wolf et al., 1998). The Large Zenith Telescope, a 6-meter 



liquid-mirror telescope, will provide photometric redshifts 
for more than 1 million of galaxies, thanks to a series of 40 
medium-band filters specially designed, and will be able to 
select quasar candidates quite efficiently. For the purposes 
of this paper, even a reduced set of broad band photomet- 
ric data on a given object shall be considered as a very 
low resolution spectrum (i2~ 1!), which could be fitted 
by synthetic templates according to the standard meth- 
ods commonly used on photometric redshift estimates. 
Assumptions on the morphology of the objects could be 
avoided in principle, in such a way that both point-like 
and spatially resolved sources could be examined through 
the same pipeline. Furthermore, it is interesting to have a 
prior idea on the redshift range of each candidate, in or- 
der to perform the spectroscopic follow up in the suitable 
wavelength range (visible or IR). 

The software described in the present paper has been 
developed for the needs of preparing the quasars/ AGN 



sample to be assembled in the future VIRMOS survey. 
It can however be applied to all multicolor catalogues. 
The forthcoming VIRMOS spectroscopic galaxy survey 
will consist in 2 surveys: a deep one and a shallow one. 
The VIRMOS-SHALLOW spectroscopic sample (statisti- 
cally 1/4 to 1/3 complete to I < 22.5, no other preselec- 
tion) will ensure random completeness, without any color 
or morphology preselection. The VIRMOS-DEEP survey 
(complete for all objects with 22.5 < I < 24) will provide 
the opportunity of assembling a quasar sample avoiding 
the drawbacks of any preselection, therefore basically free 
of the usual biases (colors, redshift etc). So far, very few 
such samples are available: three examples come from the 
CFRS (1=22.5, Schade et al., 1995), the Faint Galaxy Red- 
shift Survey (B=24, Glazebrook et al., 1995) and more 
recently Cohen et al. (1999) to R=24 in the HDF, involv- 
ing, however, very few quasars. We are examining various 
possibilities of assembling a third quasar sample in the 
limits of the VIRMOS-SHALLOW survey, trying to get 
as much as possible of the spectra of the several thou- 
sands of quasars present in this field by selecting photo- 
metric candidates. Comparing the three VIRMOS samples 
should provide a unique opportunity of examining the bi- 
ases introduced by a photometric and/or morphological 
preselection on the quasar samples. 

3. Description of the method 

Our method used here is closely based on hyperz, a pub- 
lic photometric redshift code presently under development 
at the Observatoire Midi-Pyrenees. Details on hyperz will 
be given in a forthcoming paper (Bolzonella et al., 2000; 
see also Miralles & Pello, 1998; Pello et al., 1999). It is 
basically a standard SED fitting procedure: the observed 
photometric SED of a given object, obtained through n 
filters, is compared to the SED computed for a set of tem- 
plate spectra. The aim is to find the best fit between the 
observed and the model photometry through a standard 
X 2 minimization procedure: 



X 



E 

; = 1 



{Fibs - F^ od (z)) 2 



where F* bs and i^ od are respectively the observed and 
the template fluxes in the i band, and Oi is the error on 
the observed flux in this band. Fluxes are normalized to an 
arbitrary reference filter. This is a quite general method. 
The photometric redshift is the redshift z corresponding to 
the minimum \ 2 value. As expected, the accuracy on the 
photometric redshift determination depends strongly on 
the set of filters used (number and wavelength coverage); 
this is illustrated in the next paragraph. 

In order to apply this method to quasars, the template 
set is made by series of simulated quasar spectra. These 
spectra have been constructed by varying the slope of the 
power-law spectra (spectral index) in the optical (between 
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0.0 and 1.0, in 3 steps) while keeping the ultraviolet in- 
dex, auv, constant at 1.76 (a value compatible with Wang 
et al., 1998). Emission lines are included (Ly Q , Ly^, CIII, 
CIV, Mgll, SilV, H a , Hp and H 7 ) with gaussian profiles 
and typical intensities (e.g. Peterson, 1997), as well as the 
small blue bump, centered at 3000 A. The redshift value is 
set between and 7. Totally, there are 213 quasar spectra, 
with 71 redshift steps (dz — 0.1), with 3 simulated spec- 
tra on each step, obtained by varying the optical spectral 
index. The Ly Q forest has been modelized according to 
Madau (1995), while no reddening has been included. We 
also included the stellar library of Pickles (1998), which 
contains 108 spectra of stellar types ranging from O to 
M, including dwarfs and giant stars. Additionally, we in- 
cluded 22 white dwarf and carbon stars spectra, which are 
known to contaminate the quasar candidates catalogues. 

When applying the method to a given sample, the dis- 
crimination between stars and quasars is performed in two 
steps. A first distinction between stars and quasars is made 
by identifying as star (quasars) the objects that show a 
better fit to one of the stellar (quasar) spectra. A further 
selection on the remaining sample can be made, based on 
the value of the reduced x 2 ; with stars identified as ob- 
jects excluded as quasars at the 95% or 99% confidence 
levels. In the present paper, the terms "% 2 selection" or 
"no \ 2 selection" will apply for such a distinction. When 
an object is identified as quasar candidate, a photometric 
redshift estimation is automatically provided. The num- 
ber of degrees of freedom, v, is equal to the number of 
filters minus 2: we lose one degree by normalizing the flux 
to a reference filter and one degree because of the "best 
fit" procedure. In this way, the separation between stars 
and quasars is reduced to a one-dimensional fit between 
the data (the observed SED) and the spectral templates. 

As expected, and as confirmed by the simulations, the 
quality of the quasar/star discrimination depends on the 
filter combination, as it happens for the photometric red- 
shift estimate in the case of galaxies. Not only the number 
of filters is important, but also the wavelength range cov- 
ered. Another important point is the width of the filters 
which must be suited to the width of quasar spectral fea- 
tures. A realistic determination of the photometric errors 
is the key to a successful issue for this method. 

4. Simulations 

The first test of this method has been done on mixed 
(quasars + stars) simulated catalogues. In this paragraph 
we intend to show the importance of the set of filters used 
for optimizing the results. First, let us define an arbitrary 
and useful quantity which is aimed to quantify the effi- 
ciency of this procedure. Let N c be the number of quasar 
candidates resulting from the identification technique, Nf 
the number of actual quasars found as candidates (real 
quasars effectively selected), and iV e the number of ex- 
pected quasars. In principle, N e is a quantity which can 



be estimated from previous surveys. We define the "effi- 
ciency of the method as follows: 
efficiency= completeness x confirmation rate, 
where completeness= Nf/N e and confirmation rate = 
Nf/N c . Obviously, both the completeness and the con- 
firmation rate take values in the interval [0,1] and so does 
the efficiency. 

The number of stars in our simulated catalogues was 
given by the Galaxy model of Robin et al., 1995 (R95). The 
star counts correspond to the six Deep Multicolor Survey 
(hereafter DMS) fields, since we primarily tested our re- 
sults with the results of this survey. The fields with their 
surface (in square arc minutes) and the estimated star 
numbers are given in table [l| Quasar numbers were based 
on previous surveys (Hartwick & Schade, 1990, hereafter 
HS90). Catalogues contain ~100 quasars and ~3000 stars, 
with a magnitude limit of mf,=22.2 and photometric errors 
scaling with magnitude, with typical values of the order 
of 0.1 magnitudes. Quasars are uniformly distributed in 
redshift, in order to obtain reasonable statistical results 
at all redshifts. 

Table 1. The six DMS fields, galactic coordinates, sur- 
faces and estimated star numbers. 
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Photometry in broad-band filters has been simulated, 
covering the wavelength range between ~ 3000 A and ~ 
20000 A: U, B, V, R, I, H, and K. Broad-band filters have 
been used for practical reasons, in particular because they 
directly apply to the VIRMOS survey, which will make 
use of wide-band filters only, although narrow-band filters 
would give better results (Peri, Iovino & Hickson, 1997). 
40 different filter combinations have been tested in the aim 
of selecting quasar candidates. Fig. ||a and b illustrate the 
completeness (solid line), confirmation rate (dashed line) 
and efficiency (dash-dotted line) for these 40 combina- 
tions, without any x 2 selection and with selection at 99% 
confidence level, respectively. The filter combinations in 
fig.[l]a have been classified by increasing completeness. The 
same order has been used in fig. |l]b, to demonstrate the 
effect of a x 2 selection on the completeness and the other 
parameters. Table ^| shows the filter combination which 
corresponds to the abscissa values. 

A closer look at fig. |l| and Table ^ reveals that the 
filter combination is, most of the time, more important 
than the number of filters. In general, filters like U or 
the infrared filters H and/or K are very useful: they add 
much photometric information constraining the number 
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Fig. 1. completeness (solid line), confirmation rate (dashed line) and efficiency (dash-dotted line) for 40 filter combi- 
nations, a) without any \ 2 selection and b) with selection at a 99% confidence level. 

Table 2. Filter combinations used in the fig. |l|a and |l]b, classified by increasing completeness in the case of no \ 2 
selection. 



No 


filters 


No 


filters 


No 


filters 


No 


filters 


No 


filters 


1 


VRI 


2 


VRH 


3 


BRH 


4 


UBV 


5 


BR1 


6 


BVRI 


7 


BVR 


8 


VIK 


9 


BVRIH 


10 


BRIHK 


11 


VRIHK 


12 


BRHK 


13 


BVI 


14 


BRIK 


15 


VRIK 


16 


BVK 


17 


RIHK 


18 


BVRIHK 


19 


UBR 


20 


UBH 


21 


BVRIK 


22 


UBVR 


23 


UBRI 


24 


BVIK 


25 


UBVRI 


26 


UBVRIK 


27 


UBVRIH 


28 


UBRHK 


29 


UBHK 


30 


UBVRK 


31 


UVIK 


32 


UVRIK 


33 


UBVRIHK 


34 


UVK 


35 


UBRIHK 


36 


UBVIHK 


37 


UBVRHK 


38 


UVIHK 


39 


UVRIHK 


40 


UBRIK 



of candidates. Neighboring filters like B, V, R, when used 
all together, decrease the efficiency because, as mentioned 
above, they increase the number of the degrees of freedom 
without adding any essential information. A comparison 
of the two figures reveals several trends. The completeness 
is almost not influenced by a \ 2 selection in the case of 
3-filtcr combinations (such as UBV, BRI, BVR, BVI, or 
UBR). It is, however, strongly influenced by a \ 2 selection 
when the filter U is present, in sets of 4 or more filters. 
Note that this quantity decreases significantly , when a 
X 2 selection is imposed, for filter combinations with index 
higher than 26, which are all sets containing the U-filter. 
Furthermore, the confirmation rate is generally increasing 
with increasing number of filters. The decision whether 
to make a x 2 selection or not depends on the available 
time for spectroscopy and the scientific objectives of the 
project. 

For this purpose, we made an approximative estima- 
tion of the time needed for making the spectroscopy in or- 
der to cover an area of 1 square degree, T sp , for the same 
40 filter combinations, given a detector area, S, a magni- 
tude limit, Wfc, a quasar candidate density, n q = N can d/S 
(depending on this limit as well as on the filter combina- 



tion), a number of slits, G and an exposure time, dT sp , 
with and without any x 2 selection. We calculate T sp as: 



T 



1 + int 



G 



dT s} 
S 



Fig. | illustrates T sp (in hours) for m b = 22.2, S=0.06 
square degrees, G=10 and dT sp =3 hours, without any 
preselection (solid line) and with a x 2 selection (dashed 
line). Note that the y-axis is logarithmic. As expected, 
there is an anticorrelation between T sp and the confirma- 
tion rate, as it appears from the comparison of fig. |l| and 
0. The highest the confirmation rate, the less time is re- 
quired. 

4-1- Determination of the Photometric Redshifts 

The problem of the quasar photometric rcdshift (hereafter 
z P hot) estimation is at least as complex as the one for 
galaxies. Quasar spectral features (mainly emission lines) 
are rather narrow and the overall shape of the ultra-violet 
and optical continuum is a single power low, while the fil- 
ters commonly used for photometry are wide-band filters, 
not very suitable for the quasar characteristics. All the 
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Fig. 2. Spectroscopy time (plotted on a logarithmic axis) 
per square degree versus filter combination, without any 
X 2 preselection (solid line) and with a 99% confidence level 
selection (dashed line). 



above are in the origin of the dispersed values of the esti- 
mated quasar Zphot, when compared to the spectroscopic 
ones, especially for quasars with redshifts lower than ~ 2. 
The filter combination is very important for the accurate 
determination of the z p hot- Figure |3| shows the photometric 
redshift versus the model redshift, for a catalogue of 1540 
simulated quasars, uniformly distributed over the redshift 
range [0,6], and for 4 different filter combinations. From 
top to bottom and left to right the following filter combi- 
nations are illustrated: BVRI (a), BVIK (b), UBVRI (c), 
UBVRIK (d). 

The use of the filter set BVRI is clearly a bad choice: 
these four filters cover a relatively short wavelength range 
and their effective wavelengths lie close to each other. The 
Zphot in this case is poorly determined for z < 3, while 
there is also a large dispersion at z > 5. If we replace the 
R- filter by the K-filter (figure ||b), the dispersion becomes 
more important at redshifts higher than 4, since this is 
the redshift at which Ly Q should enter the R-band, but it 
decreases the dispersion at redshifts lower than 3, due to 
the very large wavelength range covered. Adding the filter 
U to the BVRI combination decreases the dispersion at 
both high and low redshifts (figure ||c:), but one can still 
feel the absence of an infrared filter, whose importance is 
seen in figure |^d. The periodic accumulations of points at 
different redshifts are due to the passage of Ly Q from one 
filter to the other. As a general rule, the determination 
of the photometric redshifts should be rather accurate in 
the redshift range [3,5] even without the U or an infrared 
filter, but their presence is crucial for all other redshift 
ranges. 

The quality of the z p hot determination is likely to im- 
prove if real quasar spectra are used. However, this implies 
that a spectral library is required, containing some hun- 
dreds of quasar spectra, with redshifts spanning from 



to ~ 5 and with a spectral coverage spanning form the 
ultra-violet to the infrared. 

4-2. C, CR and E as a function of redshift 

The completeness, the confirmation rate and, therefore, 
the efficiency of the method depend on the redshift. The 
demonstration of this dependence is easy to show for the 
completeness, but not for the confirmation rate. The com- 
pleteness within redshift bins has the same definition as 
for the entire sample: the ratio between the number of 
real quasars with spectroscopic redshifts (hereafter z spec ) 
within the examined bin, identified as candidates even if 
their z p hot is significantly different from their z spec , over 
the estimated total number of quasars. 

Figure |J illustrates the completeness versus the red- 
shift, without x 2 selection, for 4 different filter sets: the 
solid, dotted, dashed and dashed-dotted lines correspond 
to BVIK, UBVRIH, UVK and BVRI, respectively. The 
completeness over the whole redshift range for this 4 filter 
combinations is in all cases more than 85%, as shown in 
figure ^]a. 
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Fig. 4. Completeness versus redshifts for 4 different filter 
sets: — BVIK, • • • UBVRIH, - - - UVK, BVRI. 

Figure ^ reveals that, independently of the filter com- 
bination, the completeness is always lower in the redshift 
range [2,3] due, as already mentioned, to the stellar-like 
colors of quasars. This fact could certainly increases the 
number of quasar candidates within this range but it could 
also be a possible source of incompleteness, since some 
quasars will inevitably be classified as stars. Furthermore, 
different filter combinations offer a different completeness 
at different redshifts. The filter set BVIK assures a com- 
pleteness of over 70% in the critical redshift range [2,3], 
which however remains low up to z ~ 5.5. UVK offers 
the possibility of a completeness of ^80% for z < 2, while 
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Fig. 3. Photometric redshift versus model redshift for a simulated catalogue of 1540 quasars and 4 different filter 
combinations: BVR1 (a), BVIK (b), UBVRI (c) and UBVR1K (d). 



UBVRIH results give a completeness over 90% for high-z 
quasar samples (z > 3). 

Things become more complicated if one tries to make 
an equivalent evaluation of the confirmation rate within 
redshift bins. This time, the confirmation rate is the ra- 
tio between the number of real quasars with z p h ot in the 
examined redshift bin, over the number of quasar candi- 
dates (quasars + possibly stars) with z p hot that belong to 
this particular redshift interval. Obviously, this definition 
depends on the accuracy of the determination of the z p hot 
as well as on the width of the redshift bins. Therefore, the 
dependence of the confirmation rate versus the redshift 
can be evaluated qualitatively but not quantitatively, and 
so it goes for the efficiency as well. 

If the final aim of a survey is assembling a quasar sam- 
ple with redshifts within a restrained interval, all candi- 
dates with a Zphot that lies outside this interval will, in 
principle, be rejected. In the more realistic case of an all 



redshift survey, the most interesting strategy is to sim- 
ply optimize the spectroscopy towards the optical or the 
infrared wavelengths. The spectroscopy of quasars with 
redshifts up to ~3 is much more interesting in the opti- 
cal band, while higher redshift quasars should be spec- 
troscopically observed in the infrared, since most of their 
emission lines are shifted to these wavelengths. The sim- 
ulations have shown that, although the determination of 
quasar Zphot is sometimes dispersed at z < 2.5 and z > 4.5 
(figure H) , this effect hardly translates into high-z quasars 
(z > 4) being wrongly identified as low-z ones (z < 3) or 
vice versa. 

5. Applying the method to real samples 

The method presented above has been tested on a series of 
samples found in the literature, giving rather satisfactory 
results. Depending on the sample and the available in- 
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formation, the confirmation rate and/or the completeness 
have been calculated. 

5.1. The Deep Multicolor Survey 

The Deep Multicolor Survey is an imaging survey cov- 
ering 6 fields on a total area of 0.83 square degrees on 
the northern sky at high galactic latitude, carried out at 
the 4m Mayall telescope of the Kitt Peak National Obser- 
vatory. 6 filters have been used, covering the wavelength 
range from 3000 to 10000 A: U, B, V, R\ 175 and 186, 
and a photometric catalogue of ~ 21000 stellar-like ob- 
jects has been assembled, with the aim of conducting a 
multicolor quasar search (Osmer et al. (1998) and the ref- 
erences therein). 137 stars, 49 compact narrow emission 
line galaxies (hereafter CNELGs), and 54 quasars with 
redshifts spanning from 0.3 to 4.3 have already been spec- 
troscopically identified. These objects were the 240 quasar 
candidates selected by Hall et al. according to their posi- 
tion on color-color plots, a high fraction of them (194 ob- 
jects) selected by the UVX technique (Hall et al., 1996). 
This selection favored quasars with redshifts less than 2.2. 
The redshift distribution of these objects can be seen in 
fig. ||, where the spectroscopic sample is plotted by a dot- 
ted line. 

For our purposes, we cut the catalogue at B=22.3, a 
limit at which the sample is almost 100% complete and 
only marginally contaminated (~ 8%) by galaxies, and we 
only used objects for which photometry was available in 
all 6 filters. Hence, 3 of the spectroscopically confirmed 
quasars were excluded, and we ended up with a cata- 
logue of 3720 stellar-like objects. For the magnitude limit 
and the area covered by the survey, approximately 106^22 
quasars are expected to be found among the objects of the 
catalogue (HS90), resulting a completeness of at least 80% 
for the usual multicolor selection techniques (Hall, private 
communication) and ^88% for our method. 

Table | summarizes the characteristics of the different 
catalogues examined in this section: surface covered, mag- 
nitude limits, filters used, number of quasars, and number 
of candidates for different confidence levels. The expected 
number of stars on the DMS field has been computed from 
R95. In the case of the DMS field, the number of spectro- 
scopically confirmed quasars is also given between brack- 
ets. 

Among the 3720 objects in this catalogue, we identi- 
fied 3332 as stars because their photometry shows a bet- 
ter match to the stellar templates. With our method, 6 
of the 51 spectroscopically confirmed quasars (DMS0059- 
0056, DMS1358-0054, DMS1714+5012, DMS1714+5003, 
DMS1714-4959 and DMS1358-0055) are misclassified as 
stars. Only 7 of the 97 B < 22.3 candidates, spectroscop- 
ically identified as stars, but all 34 B < 22.3 CNELGs 
belong to our candidate list, if no \ 2 selection is made. 
At a 99% confidence level, 25 CNELGs are classified as 
quasar candidates. 



To summarize, we find 388 quasar candidates for 
106 actual quasars, with an expected completeness of 
45/51=0.88 and an expected confirmation rate of (106 x 
0.88)/388=0.24, if no x 2 selection is made. The equiva- 
lent values with 99% confidence level selection are 0.88 
and 0.27, respectively, for 348 quasar candidates (see also 
Table |). 

Fig. |5| displays with a solid line the estimated redshift 
distribution of our 388 quasar candidates. The accumula- 
tion of identifications around z ~ 1 is due to the CNELGs, 
with z spec in the range [0.1,0.8] typically, identified by our 
procedure as quasar candidates with z p h ot ~ 1. The broad 
distribution at z G [2.0, 3.0] is probably due to the stellar- 
like colors of quasars at these redshifts, which makes the 
distinction between stars and quasars particularly difficult 
by photometric means, increasing the number of quasar 
candidates, but it also corresponds to the expected max- 
imum of the quasar space distribution. For comparison, 
we plot the redshift distribution of the 51 spectroscopi- 
cally confirmed quasars (dotted line). 

1 20 - 

100 - 

80 - - 

z : 1 i : 

60 - 
40 - 



0.0 0.5 1.0 1.5 2.0 2.5 3.0 3.5 4.0 4.5 
Redshift 

Fig. 5. Estimated redshift distribution of the photomet- 
rically selected candidates (solid line), spectroscopically 
determined sample of 51 DMS quasars (dotted line). 

The above results are summarized in the first 6 
columns of Table [|, where the values of the completeness 
(C) , the confirmation rate (CR) and the efficiency (E) are 
shown in columns 3 and 4. The values predicted by the 
simulations are also given for comparison in columns 5 
and 6, using the same filter combinations as for the real 
data, and keeping the filter responses as close as possible 
to those actually used. Column 2 gives the same quantities 
for the DMS and the usual multicolor selection techniques. 
Note that the values given for the completeness and the 
efficiency are the lower estimations. The two methods give 
comparable values for all 3 quantities (C, CR and E), 
with a slightly better efficiency for the multicolor tech- 
nique when a x 2 criterion is imposed. 

In order to estimate the contamination due to compact 
galaxies, we have added to the template library the spectra 
of blue star burst galaxies, obtained from the public model 
Starburst99 code by Leitherer et al. (1999), as well as a 
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Table 3. Quasars candidates for DMS, HDF Bright and Parkes selected by the z p hot method. 



Survey 


Surface 
covered 


Mag. 
limit 


Filters 


Stars 
exp 


cand 


exp 


95% 


Quasars 
Candidates 
99% 


no x 2 


conf 


DMS 


0.83sd 


B=22.3 


U,B,V,R', 


3316 ± 60 


3332 


±06122 


295 


348 


388 


51 








175,186 


(R95) 




(HS90) 


(40) 


(45) 


(45) 




HDF B 


l.OOsd 


B=21.0 


B,V,R 






30 


26 


30 


30 


30 


Parkes 




H=19.6 


B,V,R,I, 






129 


99 


107 


125 


129 



J,H,K 



Table 4. Completeness (C), confirmation rate (CR) and efficiency (E) for the samples DMS, HDF^ and Parkes and 
for the simulated catalogues, with a \ 2 selection at 99% confidence level and without any \ 2 selection. In column 2 
we give the same quantities for the DMS sample and for the usual multicolor selection techniques. 





DMS 




ot 




sim 


HDF S 


sim 


Parkes 


sim 






99% 


no x 2 


99% 


no x 2 


99% no x 2 


99% no x 2 


99% no x 2 


99% no x 2 


c 


0.80 


0.88 


0.88 


0.79 


0.86 


1.00 1.00 


0.74 0.76 


0.83 0.97 


0.52 0.81 


CR 


0.23 


0.27 


0.24 


0.20 


0.21 










E 


0.18 


0.24 


0.21 


0.16 


0.18 











delta-burst model taken from the new Bruzual & Chariot 
evolutionary code (GISSEL98, Bruzual & Chariot 1993), 
and a set of 4 empirical SEDs compiled by Coleman, Wu 
and Weedman (1980) (hereafter CWW) to represent the 
local population of galaxies (E, Sbc, Scd and Im). Spec- 
tra from the Starburst99 code were selected with solar 
metallicity and Salpeter IMF, taking different ages for an 
instantaneous burst ranging from 1 to 20 Myr. In the case 
of the GISSEL98 delta-burst, the IMF is that of Miller & 
Scalo (1979), with solar metallicity and 51 different ages 
ranging from 0.001 to 20 Gyr. CWW spectra were ex- 
tended to wavelengths A < 1400 A and A > 10000 A using 
the equivalent GISSEL98 spectra. 

When applying the method to the DMS sample of 3523 
stellar objects (we remind here that we only use objects 
with photometry in all 6 filters), using the template library 
extended with galaxy spectra, we find that ~ 80% of them 
are better fitted by stellar or quasar templates rather than 
by galaxy spectra. This fraction is only 56% for the sub- 
sample of 97 objects spectroscopically identified as stars. 
Among the subsamplc of the 51 objects spectroscopically 
identified as quasars, there are 38% that show a better fit 
to the galaxy spectra, and 80% among them correspond 
in particular to genuine young starbursts according to the 
best-fit spectra, with ages ranging between 1 and 20 Myrs. 
The result is different when the same extended template 
library is applied to the 34 B < 22.3 DMS CNELGs. In 
this case, 85% of the sample (29 objects) is well fitted 
by galaxy spectra rather than quasar or stellar templates. 
The fraction of best-fit galaxy templates corresponding to 



young starbursts is ~ 2/3. Generally speaking, a typical 
fraction between 1/4 and 1/3 of the quasar candidates is 
also compatible with young starbursts, and CNELGs are 
hardly selected as quasars, when supplementary galaxy 
templates are used. 

To summarize, adding star-bursting and early-type 
galaxy templates does not give a better solution to the 
problem of selecting quasar candidates: if star-burst can- 
didates are included in the list, the number of candidates is 
multiplied by a factor of 4 to 5. On the other side, objects 
like CNELGs will not be mixed with the quasar candidates 
but we take the risk to miss ~ 35% of the quasars, iden- 
tified as star-bursting galaxies. However, for the purposes 
of a survey, all objects will be tested through the same, 
general pipeline. Some objects (like quasars and compact 
blue galaxies) will inevitably end up as candidates in more 
than one list, depending on the different selection criteria. 

5.1.1. Quasar Photometric Redshifts 

Figure |6| left presents the z p hot versus z spec for the 51 DMS 
quasars. For 23 of them, \z p hot — z spec \ < 0.1. We note an 
important dispersion for redshifts lower than typically 2.5: 
quasars with a z spec G [1,2] are better fitted by templates 
of quasars with z S [0, 1]. This is also the tendency arising 
from simulations and shown on figure^ (upper right), for 
the filter combination UBVRI. 

Figure ^ right presents the z p hot versus z spec for a new 
sample of 46 spectroscopically identified quasars (Hall et 
al., 2000). 17 of them have \z p hot — z spec \ < 0.1. This time 
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z_spec z_spec 

Fig. 6. Photometric redshifts versus spectroscopic redshifts for the 51 quasars of the DMS sample (left) and for the 
new sample of 46 spectroscopically confirmed quasars (Hall et al., 2000) (right). 



we note that quasars with a z spec 6 [0, 1] are better fitted 
by quasars with z € [1,2], and this trend is also seen on 
the simulations (figure |3| upper right) for the same set of 
filters. 

5.2. HDFsright and Parkes 

The completeness has also been tested on two recently 
published quasar catalogues. Liu et al. (1999) published 
a list of 30 spectroscopically confirmed quasars with their 
photometry in B, V and R. These objects have 17.6 < 
B < 21.0 and 0.44 < z < 2.98, and were selected on a 
one square degree field centered on the HDF North. They 
were firstly identified as quasar candidates by their (B-V) 
(V-R) colors, along with 31 other candidates that, after 
the spectroscopy, turned out to be stars. 

With our method we identified all 30 quasars as quasar 
candidates. However, the limited number of filters did not 
allow us to reduce the number of stars selected as quasar 
candidates, 30 out of 31 were in our own list. 

We also applied our method to the Parkes subsamplc 
published by Francis et al. (1999). Only point-like objects, 
with the photometry in all 7 filters (B, V, R, I, J, H, 
K n ) were kept, thus 129 out of the 157 quasars of the 
initial sample remain in our catalogue. At a confidence 
level of 99%, 107 among these objects were classified as 
quasar candidates (completeness=0.83). If no \ 2 selection 
is applied, the completeness is equal to 0.97. 

Columns 7 to 14 in Table ^ summarize these results. 
Simulation results are also given for comparison. Note that 
in both cases, the completeness for the real samples is 
found to be even higher than expected according to the 
simulations. 



Confirmation rate and efficiency cannot be computed 
for these two samples, since the entire point-like object 
catalogues (including stars) are not available. 

5.3. ElS-wide (patch B) 

An attempt to apply our method to the ElS-wide (patch 
B) data was made. The results are approximatively the 
same as those obtained in multicolor procedures. This was 
somewhat expected: the multicolor approach proposed 
here, when applied to a 3 filter data set, is equivalent to 
a color-selection technique. Note that 3 filters have been 
used to carry out this survey: B, V and I. On a (B-V) 
(V-I) diagram, quasar evolutionary tracks up to a redshift 
of ~ 3.3 and the Main Sequence stars mix together. 



6. Discussion 

The quasar multicolor selection technique presented in 
this paper has some advantages as compared to the "tradi- 
tional" color-color methods. First, it has the advantage of 
selecting quasar candidates even at redshifts where quasar 
and stellar colors are very much alike (2.5 - 3.0). However, 
the efficiency varies with the redshift. Independently of 
the filter set, some quasars with redshifts lying in the in- 
terval [2.5,3] will be skipped. Furthermore, the number of 
quasar candidates with estimated z p hot within this range 
will be far higher in comparison to other redshift bins. The 
set of filters used plays an important role in these results. 
Adding a filter is not necessarily equivalent to adding in- 
formation and does not always improve the results, but 
imposes more stringent constraints by increasing the num- 
ber of degrees of freedom. For this reason, a x 2 selection 
does not improve the results in all cases. 
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This technique gives also an estimate of the (photo- 
metric) redshift of the candidates, and it is thus able to 
improve the spectroscopic follow up of a given sample. In 
some cases, which depend (again) on the filter combina- 
tion and the redshift of the objects, this determination can 
be accurate (\z p hot — z spec \ < 0.1). In most cases, quasars 
at redshift z < 3 will not be misidentified at z > 4 and vice 
versa, therefore the spectroscopy can be reliably targeted 
towards the optical or the infrared bands. 

For the purposes of a survey aiming in assembling 
a quasar catalogue with a minimum cost, adding star- 
bursting and early-type galaxy templates does not appear 
as a priority. On the contrary, this addition will cause 
either the increase of the candidates by a factor of al- 
most 5 (if star-burst candidates are considered as possible 
quasars) or the decrease of the completeness ~ 40% (be- 
cause an important fraction of real quasars will be bet- 
ter fitted by star-burst galaxies templates). For practical 
purposes (e.g the VIRMOS survey), when all the sources 
examined through the same pipeline, a fraction of objects 
identified as quasars will be also found within the "star- 
burst" sample, but this does not affect the present selec- 
tion. 

The same method can be extended towards different 
directions. First, it can be adapted to particular objects 
like very red quasars, by adding the suited templates. 
However, this has to be done thoroughly, to prevent from 
a possible and undesirable increase of the quasar candi- 
dates. It can also be used on spatially resolved objects, like 
Seyfert galaxies, because by construction, the software is 
independent of the morphology of the objects. This can 
be done by adding a stellar component to the AGN sim- 
ulated spectra already used. Last but not least, it can be 
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extend id to very high rcdshifts (higher than 6), after an state. edu^posmer/DMS/ 



appropriate modeling of the intergalactic absorption. 

As already mentioned, this software has been devel- 
oped for the needs of preparing the quasars/ AGN sample 
assembled in the future VIRMOS survey. However, the 
results and conclusions of this paper apply to all mul- 
ticolor surveys in general. A cross-identification between 
optical and X-ray sources can be extremely useful, when 
X-ray data are available, and can increase the efficiency, 
as quasar candidates with a compact X-ray counterpart 
will have an even higher chance of being real quasars. 
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